Finite-size scaling exponents and entanglement in the two-level BCS model 
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We analyze the finite-size properties of the two-level BCS model. Using the continuous unitary 
transformation technique, we show that nontrivial scaling exponents arise at the quantum critical 
point for various observables such as the magnetization or the spin-spin correlation functions. We 
also discuss the entanglement properties of the ground state through the concurrence which appears 
to be singular at the transition. 
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Since its experimental discovery in 1911 by Kamer- 
lingh Onnes, superconductivity has been the object of 
intensive research. More than 45 years elapsed before 
Bardeen, Cooper and SchricfFer (BCS) gave it a theo- 
retical foundation 0. The revival of interest for super- 
conductivity in the last two decades originates mainly in 
the inability of the BCS theory to explain neither high- 
Tc superconductivity nor finite-size effects in nanoscale 
grains The effect of discreteness in the energy spec- 
trum of nanograins has been studied in the reduced BCS 
model [see Eq. Q], whose exact solution was obtained 
by Richardson in 1963 Sfl 0] and whose integrability 
was proved only recently [g- 

In this Communication, we focus on the two-level re- 
duced BCS model, which displays a second-order quan- 
tum phase transition to a superconducting state, for a 
finite value of the electronic attraction. At the criti- 
cal point, we show that the spectrum and correlation 
functions possess nontrivial finite-size scaling exponents 
as already suggested for the ground-state energy 0, ■ 
Following the same line as for the Lipkin-Meshkov-Glick 
(LMG) model |l[l3, we combine a 1/A'^ expansion, the 
continuous unitary transformations (CUTs) technique, 
and a scaling argument, to exactly determine these ex- 
ponents. Our results are supported by a numerical in- 
vestigation of the finite-size effects. In a second step, we 
discuss the entanglement properties of the ground state 
via the so-called concurrence [ll| . Using a standard map- 
ping of the reduced BCS model onto a spin system w e 
show that, as in one-dimensional spin chains [l^ Il3j |. 
this concurrence displays some singular behavior at the 
transition point. In the present case, one may, however, 
distinguish between two cases depending on which level 
the two spins considered for the concurrence belong to. 

We consider the reduced BCS Hamiltonian 
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where c] ^ and c. 



'i,± 



are fermionic creation and annihi- 



lation operators obeying the anticommutation relation 
[Cj g., c|, ^,]+ = Si^ir6a,a', 5 is a positive coupling constant 



that we set to unity, N is the total number of states 
{j = 1, . . . , N) that is assumed to be even, and the 1/A^ 
factor ensures that the thermodynamical limit is well be- 
haved. For the sake of simplicity, wc furthermore restrict 
ourselves to the subspace where all states j have occupa- 
tion number or 2. This subspace is not coupled to its 
complement by the BCS Hamiltonian. Via the mapping 
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of fermionic to spin operators, the two-level Hamiltonian 
at half-filling becomes 



H = -h{Si, 
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We have denoted ±/i the energies of the two levels 
(/i > 0), which are interpreted as magnetic fields in the 
spin language. We have introduced total spin operators 
Si(2)a = Z]»e/i(2) ^a/^' ^1(2) being the set of states i with 
energy +h {—h), and Sa = Sia + S2a (where a — x,y,z). 

Note that we have also used the conservation of the 
number of fermions, which reads 5"^ = in the spin lan- 
guage. The Hamiltonian Q further commutes with the 
total spin operators and S^. Here, we focus on the 
maximum spin sector to which the low-energy states be- 
long to. 

In the spin language, the two-level BCS Hamiltonian 
thus describes the physics of two coupled XX models 
with infinite range (constant) interactions, and embed- 
ded in transverse magnetic fields of opposite directions. 
In the thermodynamical limit, the Hamiltonian JSjl un- 
dergoes the well-known BCS mean-field quantum phase 
transition at /i = 1, that we now briefly describe. Re- 
placing the large spin operators by their classical value 
(fc = l,2) 

Sfc = ^(sin6'fcCOS(/)fc,sin6'fcSin(?!)fc,cos6'fc), (4) 
and minimizing Q with respect to 9 and 0, yields the 



solution 



— 0,02 — vr) in the symmetric phase h > 1, 
and {91 = arccos ft,, 6*2 = tt — 91) in the broken phase h < 
1. The spontaneous breaking of the rotational symmetry 
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TABLE I: Scaling exponents for the ground-state energy per 
spin eo, the gap A, the magnetization (Siz), and the two- 
point correlation function (Sl^) and {SixS2x)- 



along z for h < 1 gives rise to infinitely degenerate ground 
states {(pl = (j)2 can take any value). 

We compute the finite-size scalin g ex ponents following 
the same strategy as in Refs. |^ andll5 We consider the 
symmetric phase h > 1, and use the bosonic Holstein- 
Primakoff representation 0| of both spins around the 
mean-field ground state 
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S2+ = iS2-y = iN/2y/%^ {l-2b^b/N) 
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One then inserts these in Eq. |(2Jl, and expands the square 
roots to the lowest order needed to compute the expo- 
nents, i.e., (l/iV)^. The Hamiltonian is not diagonal and 
contains terms creating or destroying one a and one b 
excitations, for example a^6^ -I- ab. Since it is quartic, 
it cannot be simply diagonahzed by a Bogoliubov trans- 
form, but one can use the CUTsIm lia lia to perform 
this task (see also Refs. 18 andllO|for details). Starting 
from the initial Hamiltonian H{Q) = H, one considers a 
unitary equivalent Hamiltonian H(l) satisfying the flow 
equation 



d,H{l) = m,H{l)]. 



(8) 



Tjil) is the generator of the unitary transformation, cho- 
sen to make the final Hamiltonian H{1 = oo) diagonal. 
For the problem at hand, the simplest generator is the 
so-called particle-conserving generator [ijj, with parti- 
cle number operator Q = ria + ni,, where Ua = a^a and 
nb — b^b. With this choice, the final Hamiltonian is poly- 
nomial in Ua and n;,. To compute correlation functions, 
one has to follow the flow of spin operators. Sz = Ub — Ua 
is found to have no flow, so that eigenstates of the final 
Hamiltonian must satisfy Ua — rib to fulfill Sz ^ 0. The 
ground state (first excited state) of H{1 = oo) is the state 
with Ua = 0(1) boson. 

Concerning the spectrum, we focus on the ground-state 
energy per site eo and the gap A. For the magneti- 
zation, one has (S'2z) — ~(5'iz), and the rotation in- 
variance around the z axis implies (Skx) — (Sky) = 
for k = 1,2. Finally, all spin-spin correlation functions 
{SjaSkfS + SkpSja) with j, k e {1, 2} and a, /? e {x, y, z} 



either vanish or can be deduced from (Sf^) and (<Sj^S'2x)- 
As in our earlier works on the LMG model |9|, [lOj , the 
flow equations can be integrated exactly, and the five 
quantities of interest are found to behave as 
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where E{h) — h{h — 1). 

Let us now outline the argument already used in Refs.0 
andliolto compute finite-size scaling exponents. The 1/N 
expansion of any physical quantity <I> considered here has 
a simple structure. It consists of two contributions which 
are, respectively, regular (reg) and singular (sing) when 
h approaches the critical point. Schematically, one has 



^M{h) = <i>^«(/i) + <^7Hh). 



(14) 



Nevertheless, it is clear that no divergence can occur at 
finite N for these quantites or its derivatives with respect 
to h, even at the critical point. This basic fact straight- 
forwardly leads to the scaling exponents. Indeed, a close 
analysis of the singular part <I>™^(/i) allows us, in the 
vicinity of he, to write it as follows: 



S(/i)«- 



TV"* 
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where JF$ is a function that only depends on the scaling 
variable A^S(/i)3/2. To compensate the singularity com- 
ing from S(/i)^*, one thus must have JF$(x) ~ 2;-2C4./3 
so that «'57^(1) - Ar-(»*+2?*/3)^ The scaling exponents 
corresponding to l|^- H13|l are listed in Table HJ 

We emphasize that here, we have approached the crit- 
ical point from the symmetric phase (h — > 1"^). How- 
ever, we could, as for the LMG model have reached 
it from the broken phase. As pointed out by Richard- 
son [i^l who derived the exact solution in this regime, 
the 1/N developments of the ground-state energy and 
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FIG. 1: Behavior of the ground-state energy per spin eo and 
the gap A as a function of the system size TV (logj — logj plot) 
at the critical point compared to analytical results. 



FIG. 2: Finite-size scaling of the magnetization and the cor- 
relation functions (log2 — logj plot) at the critical point. A 
clear power-law behavior is observed with a nontrivial scaling 
exponent 2/3 predicted by the present method. 



the occupation number (which is essentially the magne- 
tization in the spin language) are also singular at the 
critical point. We have checked that the scaling argu- 
ment given above also predicts eQ"^(iV) ~ N~^/'^ and 
4(S'iz)^'"s/7v ~ Ar-2/3 from this side of this transition, 
i.e., when ft, — s- 1~. Recently, the nontrivial scaling expo- 
nent of eo has been observed numerically and analyt- 
ically derived using a spin coherent state representation 
[3 but, to our knowledge, the other scaling exponents 
given m Tabled have never been discussed. The present 
results can be compared to those recently obtained in the 
LMG model HII^I- For the isotropic LMG model, which 
has the same interaction term (5*^ + Sy) as the Hamil- 
tonian (O, scaling exponents are trivial since the 
expansion is not singular at the critical point. However, 
for the anisotropic LMG model, similar exponents are 
found (multiple of 1/3) except that x and y directions 
have different exponent. 

To check the validity of the present approach, we have 
performed numerical diagonalizations of critical finite- 
size systems, with up to N = 2^^ spins for eo and A 
and up to A'^ = 2^'^ spins for the magnetization and the 
correlation functions, which require the knowledge of the 
eigenstates. The results are depicted in Figs. ^ and 13 
and show an excellent agreement with the analytical pre- 
dictions. Note that the regular part $^^(1) has been 
substracted to underline the nontrivial scaling behavior. 

Let us now discuss the entanglement properties of the 
ground state. Here we focus on the concurrence [ll| . 
which characterizes the entanglement between two spins. 
As in the LMG model , in the thermodynamical limit 
the ground state becomes a completely separable state. 
Thus the nontrivial properties of the concurrence are en- 
coded in the finite N corrections and one has to consider 
the rescaled concurrence Cr = {N — 1)C. For the two- 
level reduced BCS model, one has to distinguish between 
two cases: (i) both spins belong to the same subset (/i 
or 12); {a) each spin belongs to distinct subsets. In both 



situations, one can express the rescaled concurrence as a 
function of the observables previously calculated. Gen- 
eralizing for our purpose the results given in Ref. one 
gets 



Cp^ — Cp^ = 2max(0,w 
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The superscript refers to the subset to which the spins 
considered belong to, and 
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Using expressions (lllll - (|13|) truncated to order (1/A^)^, 
one can compute the thermodynamical limit of the 
rescaled concurrence in the symmetric phase (h > 1), 
and one finds 
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In the broken phase (0 < /i < 1), one has to use the 
Holstein-Primakoff representation around the mean-field 
ground state. We found out that it is then enough to 
perform a Bogoliubov transform to compute the concur- 
rence. Indeed, although such a simple calculation fails 
to provide the full (1/A^)^ contribution to the spin ex- 
pectation values (as explained in Ref. for the LMG 
model), the unknown parts of the (1/A^)^ contributions 
cancel each other when computing the thermodynamical 
limit of the rescaled concurrence, and one obtains 



Cr ' = max (0,2- ^ ^ 
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FIG. 3: Rescaled concurrence of the ground state as a function 
of the magnetic field for A'^ = 32, 64, 128, 256, and oo. In the 
thermodynamical limit a cusplike singularity appears at the 
critical point. 

The rescaled concurrence is depicted in Fig. Ofor several 
system sizes as well as in the thermodynamical limit. Let 
us also mention that the finite-size scaling exponent for 
can be computed using the same scaling argument 
as previously and equals 1/3 instead of 2/3 as expected 
from the scaling of the observables. 

As in the LMG or the Dicke model '23'], it is interest- 
ing to note that displays a cusplikc behavior at the 
critical point, whereas it is a smooth hmction for h ^ \. 
More interestingly, the max function used in definition 
(PI leads to Cr'^ = for /i > ^/3/2. ft is worth noting 
that this value of the magnetic field does not play any 
special role in the phase diagram, whereas it naturally 
arises in the entanglement analysis. However, if we do not 
consider the max function, C^'^ diverges at the critical 
point. In the zero-field limit, the ground state is simply 
the Dicke state corresponding to 5* = A^/2 and 5*2 = 0, 
whose rescaled concurrence equals 1 [13, 0|- Note that, 
in this limit, the distinction between subsets 1\ and I2 
becomes irrelevant so that both rescaled concurrences are 
identical. In the infinite h limit, the ground state is the 
separable state IV^o) = (®ie/ilT)i) ® (®je/2li)j)7 which 
has also a vanishing total magnetization (Sz — 0) but 
which is not an eigenstate of S^. However, it is clear that 
the concurrence of such a state is exactly zero whatever 
the two considered spins. 

The method used in the present work is certainly 
well suited to tackle many similar problems for which a 
semiclassical description of the thermodynamical limit 
is exact as illustrated here or in the LMG model. An 
interesting issue is to understand the key ingredients 
that allow for exact solutions of the flow equations. In 
particular, it would be of special interest to analyze, 
along the same line, models where some regions of the 
phase diagram are known to be chaotic and others 
integrable. Finally, let us note that our method may also 
be useful to compute the von Neumann entropy which 
displays some interesting features in quantum critical 
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